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Abstract 



We study an extension of the gravity dual to a perfect fluid model found by Janik and 
Peschanski. By relaxing one of the constraints, namely invariance under reflection in the 
longitudinal direction, we introduce a metric ansatz which includes off-diagonal terms. We 
also include an i?-charge following Bak and Janik. We solve the Maxwell-Einstein equations 
and through holographic renormalization, we show that the off-diagonal components of the 
bulk metric give rise to heat conduction in the corresponding CFT on the boundary. 
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1 Introduction 



It is intriguing that strongly interacting particles produced in heavy ion collisions exhibit perfect 
fluid characteristics which are appropriately described by a hydrodynamic model. The hydrody- 
namic behavior in the context of Quantum Chromodynamics (QCD) was studied by Bjorken in 
the early eighties [1]. It remains a challenge to derive such behavior from first principles in QCD. 
On the other hand, recent experimental results at RHIC have provided experimental evidence of 
a hydrodynamic form of strongly interacting matter [2]. 

Much has been learned about gauge theories in the strong coupling regime through the AdS / CFT 
correspondence [3] [4] which enables us to understand these gauge theories by studying their grav- 
ity duals in AdS space with one additional dimension. Strictly speaking, this correspondence 
applies to the maximally symmetric M = 4 super Yang-Mills theory. Nevertheless, its applica- 
tion to gauge theories with less symmetry, including QCD, has been considered. In particular, 
regarding heavy ion collisions, a gravity dual has been used to extract information about jet 
quenching [5], transport coefficients [6], the fireball produced at RHIC [7], as well as a bound on 
rj/s[S}. 

In an interesting work, Janik and Peschanski [9] discussed a solution to the Einstein equations 
in the bulk AdS space in the limit r — > oo, where r is the proper time in the longitudinal plane. 
By demanding regularity in the bulk, they showed that the acceptable solution corresponds 
through holographic renormalization [12] to a perfect fluid on the boundary of AdS. The work 
was furthered by Bak and Janik [10] who studied the Maxwell-Einstein equations in the bulk 
with a conserved R-charge. 

Here we extend the results of [9] [10] by relaxing the constraint of invariance under reflection in 
the longitudinal direction. This allows us to include off-diagonal terms in the bulk metric ansatz. 
We solve the Maxwell-Einstein equations in the bulk in the limit r — > oo. We obtain an exact 
solution and show, through holographic renormalization, that it corresponds to a non-viscous 
fluid of non-vanishing chemical potential. The novel characteristic is a temperature gradient in 
the longitudinal direction. We study the thermodynamic properties of the fluid and their relation 
to the form of the bulk metric. 



2 R- Charged Perfect Fluid 

We are interested in understanding the behavior of a fluid described by a gauge theory in a 
four-dimensional space spanned by coordinates x M , (/i = 0,1,2,3). Following [1], we introduce 
the proper time r and rapidity y on the longitudinal plane, defined by 

x° = TCOshy , x 1 =rsmhy (1) 

The transverse coordinates will be denoted by x 1 - = (x 2 ,x 3 ). The gravity dual of the four- 
dimensional theory will be five- dimensional. Let z denote the fifth dimension. We shall solve the 
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Maxwell-Einstein equations in the bulk 

Rap + ^ + ^F^F^ - = (2) 

V a F Q " - -^ e ^ 5 F aP F l5 = (3) 
48^ — g 

Following [10], we adopt the metric and field strength tensor ansatz 

ds 2 = ^[-e a ^dr 2 + r 2 e b ^dy 2 + dx L2 + e d ^dz 2 ] (4) 
F ZT = K(t,z) (5) 

which is the most general bulk metric obeying boost invariance, symmetry under reflection in the 
longitudinal direction (y — > —y), plus translational and rotational invariance. It is convenient to 
introduce the coordinate 

v = — (6) 

in terms of which the metric (4) reads 

ds 2 = -JL [_( e «(^) - e d ^r] 2 v 2 T 2 ^)dT 2 + r 2 e b ^dy 2 + dx ±2 + e d ^r 2r >dv 2 ] 

+ 2^)*^ (7) 
F w = F(r,u) (8) 

Substituting this ansatz into the Maxwell-Einstein equations, we obtain the leading behavior in 
the r -> oo limit [9] [10] 

, = 1/3, F(r,v) = ^ 

a{y) = In (\ + av 4 + |^ , 6(u) = 0, = - In (\ + av 4 + |^ 6 ^ (9) 

in terms of arbitrary parameters a and q. 

The above bulk metric may be related to the vacuum expectation value of the stress-energy tensor 
of the corresponding gauge theory on the boundary through holographic renormalization [12]. 
This is done as follows. The metric (4) is transformed to that of the general asymptotically AdS 
metric in Feffer man- Graham coordinates 

ds 2 = 9 ^ dX J + ^ (10) 
z 2 

Near the boundary at z = we may expand 

9r = $ + * l 9]2 + *'$ + ■■■ (11) 
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where gf} = r]^ and g$ = 0. The next coefficient, g$, is proportional to the vacuum expecta- 
tion value of the stress-energy tensor. Thus, 



N 2 

{Tuv) = ^} (12) 



It is can be seen from the form of g$, for the special value (9) of r], that the stress-energy tensor 
corresponds to that of a perfect fluid 

T af s = (e + p)u a Ufs + pr] a/3 (13) 

obeying the equation of state p — |e (tracelessness due to conformal invariance). Further, 
from the explicit form of the metric it follows that the energy density and temperature fall off, 
respectively, as 

T = T H = — (l- q ^Pj r-V3 (15) 

7TV h V 24 J 

where Vh is the position of the horizon (obeying 1 + av\ + q 2 v^/12 = 0). 

After performing a detailed thermodynamic analysis one finds the entropy density (s), charge 
density (p), and chemical potential (/x), respectively, 

Swi A'o 2 D" N?q qvl 

S= V^r^ P= VL; = 8A' M = >M*)-Ar<* = 0) = ^ (16) 

Expressing the energy density in terms of s and p, 

/ \ 3s 4 / 3 / 4vr 2 p 2 \ 

one verifies that the chemical potential is conjugate to p (p = de/dp). 

We must take note that that the geometry is evolving and exact notions of temperature and 
entropy are ill defined. However, it is assumed that the approximate notions still prevail. A 
discussion of dynamical horizons is given in [11]. 



3 An extension of the bulk metric 



Next, we consider an extension of the bulk metric (4) which yielded the perfect fluid model 
discussed above. To this end, we relax one of the conditions which led to the general form of 
the metric, namely the requirement of invariance under reflection in the longitudinal direction, 
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y — > —y. We may then add appropriate off-diagonal terms in the bulk metric (7) that will lead 
to new non-singular solutions to the Maxwell-Einstein equations (2). We therefore consider the 
metric and gauge field ansatz 

ds 2 = — [-(e a ^ - e d ^r] 2 v 2 T 2 ^)dT 2 + r 2 e b ^dy 2 + dx ±2 + e d ^r 2r >dv 2 ] 
+ 2 V e d ^— + 2v- 2 T X - 2r >h ll (v)dTdy 

VT 

A T = r^A (v) 

A y = r^A 1 (v) (18) 

where the second line of the metric consists of the off-diagonal terms added to the "perfect 
fluid" ansatz (7). The exponents — 1 < A — 2i],t; ,£ il < of r, where bounds are placed so that 
the associated functions do not singularly dominate, and the functions, a(v), b(v), d(v) and 
longitudinal coupling h\\(v) are to be determined by the Maxwell-Einstein equations (2). 

Substituting the modified ansatz (18) into (2), we obtain coupled differential equations for the 
various functions parametrizing the ansatz. The task of extracting the leading contribution is 
complicated by the fact that we now have four parameters determining the order of the expansion, 
rj (which is present in the perfect fluid case) as well as A, £ > and £i- For consistency of the 
expansion, we need 

< 77 < 2/5 (19) 

to be compared with the acceptable range < rj < 1 in the perfect fluid case. The restricted 
range still includes the special value (9). 

Let us first consider the v and r components of the Maxwell equations. They may be factored 
in powers of r as 

e a(v)+2b(v) v 4 A ^ v) ^ Q + 1 _ 2r) ) + ( r 2(A-l) > r A-K 1 -2-& > r 3A +?1 -4- ?0) = Q 

e o(t,)+26( w ) u3 [A' ( v )(2 + va'(v) - vb'(v) + vd'(v)) - 2vA%(v)] 

+O(r 2 ( A - 1 ),r A+ « 1 - 2 -«°,r 3A+ « 1 - 4 -« ) = (20) 

The v component is satisfied by the choice of parameter 

£o = 2 V - 1 (21) 

The r component is solved by 

A' (v) = qv e h(°W-t>(v)+d(v)) ( 22 ) 
where q is an arbitrary integration constant. 

Turning attention to the y component of the Maxwell equations, we observe that the leading be- 
havior is not unambiguously determined without further information on the parameters. Indeed, 
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in addition to the manifestly O(t ) term, there are two other terms in which the exponent of r 
is not necessarily negative. The y component reads 



2a(v)+b(v) 3 



1 ^'v 3 [A[(v) (-2 + va'(v) - vb'(v) - vd'{v)) + 2vA'[{v)] 

_ T X-Zi-l/3 e a(v)+b(v) v 3 

x [-2vA^(y)h\\(v) + h\\(v) (A' (v) (2 + ua» + vb'(y) + ud'(u)) - 2vA£(v))] 
_ r 3A-6 -7/3^3^3^) [A[ ) (t;)(2 + ^ , (t;)) -2^'^)] + C(r 2(A - 1) ,r 2 ^- 1 )) = (23) 

A similar ambiguity occurs in the and i>y components of the Einstein equations. They read, 
respectively, 

Cl ( v yn-2 + C2 ^ v ) r ^-x-i + ai (u)^||(u) + ^{v)h'^v) + /^(u) + 0(r 2 ^ +A - 2 )) = 

boivy^-^ 1 + b^)/^) + \>2(v)h\\(v) + h'^v) + 0{r 2 ^) = (24) 

where 

„.2 p —a ii A' 

Cl = ^-/^H^ 2 , c 2 = „X4 , b = ^ ha; - Ci^i) 



8(i -e d ) w yn+i+A) , y, , 

^ 1 h a , bi = h — (a 

ir v vrj 2 



a 2 = — -I( a ' + 6' + d') , b 2 = ?7 1 A -^(a / + 6 / + d / ) (25) 

v 2 vrj 2 

For consistency of the two equations (24), we must have a 2 = b 2 , which leads to the constraint 
on the parameters 

A = Ar] - 1 (26) 

The condition (26) is necessary for consistency but not sufficient. We need to ensure ai = bi, 
as well. To this end, we turn to the diagonal components of the Einstein equations which will 
determine the parameters rj and £i as well as the functions a(v), b{v) and d(v). As of now there 
is enough information to determine i] by demanding finiteness of the Maxwell scalar F 2 . We 
have 

F 2 = F af3 F a ? = -2gVe-^V 6 ^ 2 + 0(r 8r >^-\ r 2 ^^) (27) 
which forces rj (and therefore A on account of eq. (26)) to take the special value 

v = A = 1/3 (28) 

The y component of the Maxwell equations may now be used to determine the remaining pa- 
rameter £i. Since £i < 0, the exponents of r involving £i will be positive unless 

6 = (29) 

With this choice, all terms shown explicitly in the y component of the Maxwell equations are 
O(t ) (leading). 
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The rest of the Einstein equations are found to 0(t 4 / 3 ) as 

lQ(e d - 1) + \e~ a v A A' 2 + 2v (4a' + b' - a") + v 2 (-a' 2 - a'b' + a'd' - 2a") = 

16(e d - 1) - \e- a v A A% + 2*; (a' + 46' - a") + v 2 (-a'b' - b' 2 + b'd' - 2b") = 

8(e d - 1) - \e~ a v A A' 2 + v(a' + b'- d') = 
3 

16(e d - 1) - \e~ a v A A% + 2t; (a' + b' - 4a") + v 2 (-a' 2 - b' 2 + a'd' + b'd' - 2a" - 2b") = 

16(e d - 1) + \e~ a v A A' 2 + 2 (4a' - 6' - a") + t; 2 (-a' 2 - a'6' + a'd' - 2a") = (30) 

To leading order, these are identical to the perfect fluid case and are solved by (9). The off- 
diagonal component h\\(v) of the metric enters in the next-to-leading order. Once again, we note 
that the expansion is consistent for the restricted ranges of 77, £0, £1 and A. 

Turning back to the equations (23) and (24), we observe that for the functions (9), the two 
Einstein equations coalesce and combined with the y Maxwell equation Ai(v), hu(v) are uniquely 
determined at leading order. The resulting differential equations have regular solutions found 
after some straightforward algebra as 

h\\(v)=A^ + ^v^ (31) 

and 

A'Jv) = - q —v (32) 
a 

where A is an arbitrary constant characterizing the departure from the symmetric state (under 
reflection in the longitudinal direction) discussed in [10]. 



4 Hydrodynamics 

In order to understand the dynamics of the corresponding gauge theory on the boundary we 
invoke holographic renormalization [12]. To do this we must pass back to the Fefferman-Graham 
coordinates (10) and obtain the fourth order term in the expansion. Following the same procedure 
as in [10], we redefine 

* = ^(! + ^f + -) (33) 

which effectively flattens the z coordinate and from eq. (12) we may determine the vacuum 
expectation value of the gauge theory stress-energy tensor. The diagonal components maintain 
the same form that led to the previous energy density and pressure (14). The off-diagonal 
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component of the metric demands the vacuum expectation value of the gauge theory stress- 
energy tensor develop an off-diagonal component which behaves as 

N 2 A 

= ^ 7 (34) 



In order to understand how our solution relates to the gauge theory fluid, let us choose a stress- 
energy tensor which includes an arbitrary energy flux in the longitudinal direction, start by 
working with an arbitrary stress tensor with diagonal and ry components, 



/ rpTT 


rp T y 








\ 


rpry 


















rp22 







V o 








2^33 


/ 



(35) 



Recalling our definition (1) of coordinates r and y, the metric on the Minkowski space of the 
fluid reads 

ds 2 A = -dr 2 + r 2 dy 2 + (dx^ 2 (36) 

From the local conservation law 

V a T aP = d a T aP + T« X T XI3 + x T aX = (37) 

and using the Christoffel symbols — ^ — T^. and Ty y = r, we derive relations between the 
components of the stress tensor. 



Choosing (3 = r, we obtain 

d T T TT + d y T Ty + -T TT + tT vv = (38) 

T 

Setting f3 = y, we deduce 

d T T Ty + d y T yy + -T Ty = (39) 

T 

One more relation is a consequence of conformal invariance. Demanding tracelessness, we obtain 

_ T rr + T 2 T yy + + ySS = q ^q^ 

In order to match with the expected form of the stress-energy tensor from holographic renor- 
malization, we observe that the components of the stress-energy tensor to the order we are 
considering should not depend on the rapidity y or the transverse coordinates ark We may then 
immediately solve for the energy flux component T Ty , obtaining 

T- = ^ , T Ty = ^ (41) 

with C being an arbitrary constant. This behavior matches the prediction (34) of the gravity 
dual for T Ty . 
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Solving for the diagonal components T yy ,T n (i = 2,3), we obtain the stress-energy tensor in 
terms of e — T TT (energy density) and the arbitrary constant C (determining the energy flux), 



( e £ 


\ 



T 



\ 



e + \rd T e 

e + \rd T e ) 



(42) 



The functional form of the energy density e is obtained by demanding isotropy, r 2 T yy = T 22 = 
T 33 . We deduce 

£ - ^3 ( 43 ) 

matching the perfect fluid behavior (14). The equation of state is also unchanged, p = |e. For 
the energy flux, we obtain the velocity component 



u y ~ r 



rpry 



1 



r 2/3 



(44) 



Summarizing, we have obtained the following behavior of the components of the stress-energy 
tensor, 

T"=* . T, 5 



r 4/3 ' r 4/3 

T OT _ ^ T _ B 2/3 

- 3r io/3 > 1 w- 3 r 

T " = 3^ ' Tii = 3^ ^ = 2 ' ^ 

Tv = ^ , T r , = ^ (45) 

This behavior exactly matches that expected from the gravity dual. The constants B and C are 
related to the parameters of the bulk metric by 

B=™L(l+W\ C = ^A (46) 



To gain further insight into the nature of this fluid, consider the general case of a dissipative 
relativistic fluid with stress-energy tensor and current respectively given by [13] 

T afi = (p + p)u a u p + pr] aP + t afi , J a = pu a + f (47) 

generalizing the perfect fluid expression (13). p is the particle density (for a perfect fluid, J a = 
pu a ) and the dissipative pieces satisfy 

u a t a/3 = u a f = (48) 
Standard thermodynamic arguments using the relations 

e + p = Ts + pp , de — Tds + pdp , dp = sdT + pdp (49) 
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lead to the general expressions 

f0 = -r)[V a U P + u a v?Vy + (a <-> 0)\ ~ (c ~ \vj [r] af3 + uV]V 7 u 7 (50) 

j a = -*(d a + u a u%)^ (51) 

in terms of the coefficient of thermal conductivity x and shear and bulk viscosities, r\ and (, 
respectively. Assuming no particle transport, J = , we deduce 

pu = x(V + tiu%)£ (52) 

where the vectors are three-dimensional. This can be solved as an expansion in the derivatives 
of p/T, 

u = -V- + . . . (53) 

This expansion is justified because local velocities are small (\u\ <C 1). To match the behavior 
(45) of the stress-energy tensor, we assume no viscosity, thus setting 

V = C = (54) 

Furthermore, the energy flow ought to be in the longitudinal direction, thus 

u L = (55) 

We may determine the energy density and pressure (14) directly from the form of the bulk metric. 
The results are in agreement with the symmetric state considered in [10] (see section 2). The 
other thermodynamic quantities (temperature, entropy, charge density and chemical potential) 
acquire dependence on the longitudinal coordinate (as well as r) due to the heat flux. Thus, we 
obtain corrections to the symmetric case (eqs. (14) - (16)). Expanding in the y-coordinate, we 
may write 

T = T + yT 1 + . . . , s = s + ys 1 + . . . , p = p + ypi + . . . , p = p + W\ + ■ ■ ■ (56) 

where the zeroth-order contributions are in agreement with their counterparts in the symmetric 
case. Higher-order corrections may be determined from the thermodynamic relations and the 
form of the stress-energy tensor. 

where T is the ideal system's temperature and 7\ is the first order correction. The remaining 
thermodynamic terms are similarly expanded. 



Making use of (49) we obtain 



rji S T 

Pi = -Ti— , p 1 = -si — (57) 



Po Po 

Combining with (53) we find 

y X //'I ;■ t'<>\ X ( ■% Po . 



From (47), the energy flux in the longitudinal direction is 

T ry = {e + p)u T u y (59) 

To leading order, this is given by 

where we used the equation of state p = |e. 

Upon comparison with eqs. (14) - (16), (34), we deduce the temperature and chemical potential 
gradients, 

1 327T 4 X (l + g2 v 6/ 12 )2 T ' A*l (1 + ^2^6/^)2 ^ 

the entropy and charge density gradients may be found from the thermodynamic relation (17), 

. = <M^_ 40: - <MI^£ _-5/3 

1 64tt 4 x (1 + g 2 u«/12) 2 (l + g 2 u£/24) 

9 .4iV c 4 ^(1 - ^g/24) 3 5/3 , . 

^ 327r 5 x(l + g 2 ^/12) 2 (l + g 2 ^/24) 1 J 



Finally, the energy flux may be written in terms of the temperature gradient (since the pressure 
gradient is of higher order) as 

T Ty = -K T d y T (63) 

where kt is the standard definition of thermal conductivity in non-relativistic mechanics. We 
obtain [6] 

Kt = K (^£) (64) 



[rpT 

To leading order, we have 

_ 167T 2 (l + gM/12) 2 

and using the expression (61) for the temperature gradient, we easily see that the energy flux 
(63) matches (34). 



5 Conclusion 



We solved the Maxwell-Einstein equations in AdS*, for long longitudinal proper time using a 
metric ansatz which was a generalization of the proposal of ref. [10]. By relaxing the requirement 
of invariance under reflection in the longitudinal direction, we were able to add off-diagonal 
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terms to the metric. We found an explicit solution by keeping leading terms in the Einstein 
equations. Through holographic renormalization, we showed that the bulk metric corresponded 
to a non-viscous gauge theory fluid with energy flux in the longitudinal direction. We studied its 
thermodynamic properties and calculated the standard coefficient of thermal conductivity and 
the temperature, chemical potential, entropy, and charge density gradients. 

It would be interesting to relax the assumption of rapidity invariance of the bulk metric and con- 
sider more general forms of the metric by including dependence on transverse spatial coordinates 
as well. That would lead to viscous behavior as well as a more complicated energy flow. Work 
in this direction is in progress. 



Acknowledgments 

The work of J. A. and G. S. was supported in part by the Department of Energy under grant 
DE-FG05-91ER40627. C. M. was supported in part by a Faculty Professional Travel Grant from 
Rhodes College. C. M. gratefully acknowledges the hospitality of the Department of Physics and 
Astronomy at the University of Tennessee where part of this work was done. 



References 

[1] J. D. Bjorken, "Highly Relativistic Nucleus-Nucleus Collisions: The Central Rapidity Re- 
gion," Phys. Rev. D 27, 140 (1983). 

[2] P. F. Kolb and U. W. Heinz, "Hydrodynamic description of ultrarelativistic heavy-ion col- 
lisions," arXiv:nucl-th/0305084. 

[3] J. M. Maldacena, "The large-A" limit of superconformal field theories and super gravity," 
Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 (1999)] [arXiv:hep- 
th/9711200]; 

S. S. Gubser, I. R. Klebanov and A. M. Polyakov, "Gauge theory correlators from non- 
critical string theory," Phys. Lett. B 428, 105 (1998) [arXiv:hep-th/9802109]; 
E. Witten, "Anti-de Sitter space and holography," Adv. Theor. Math. Phys. 2, 253 (1998) 
[arXiv:hep-th/9802150]. 

[4] O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri and Y. Oz, "Large- N field theories, 
string theory and gravity," Phys. Rept. 323, 183 (2000) [arXiv:hep-th/9905111]. 

[5] S. J. Sin and I. Zahed, "Holography of radiation and jet quenching," Phys. Lett. B 608, 265 
(2005) [arXiv:hep-th/0407215]. 

[6] D. T. Son and A. O. Starinets, "Hydrodynamics of R-charged black holes," JHEP 0603, 
052 (2006) [arXiv:hep-th/0601157]. 



11 



[7] E. Shuryak, S. J. Sin and I. Zahed, "A gravity dual of RHIC collisions," arXiv:hep- 
th/0511199. 

[8] G. Policastro, D. T. Son and A. O. Starinets, "The shear viscosity of strongly coupled 
N = 4 supersymmetric Yang-Mills plasma," Phys. Rev. Lett. 87, 081601 (2001) [arXiv:hep- 
th/0104066]; 

G. Policastro, D. T. Son and A. O. Starinets, "From AdS/CFT correspondence to hydro- 
dynamics," JHEP 0209, 043 (2002) [arXiv:hep-th/0205052]; 

P. Kovtun, D. T. Son and A. O. Starinets, "Viscosity in strongly interacting quantum 
field theories from black hole physics," Phys. Rev. Lett. 94, 111601 (2005) [arXiv:hep- 
th/0405231]; 

A. Buchel and J. T. Liu, "Universality of the shear viscosity in supergravity," Phys. Rev. 
Lett. 93, 090602 (2004) [arXiv:hep-th/0311175]. 

[9] R. A. Janik and R. Peschanski, "Asymptotic perfect fluid dynamics as a consequence of 
AdS/CFT," Phys. Rev. D 73, 045013 (2006) [arXiv:hep-th/0512162]. 

[10] D. Bak and R. A. Janik "From static to evolving geometries-R-charged hydrodynamics from 
supergravity," Phys. Lett. B 645, 303 (2007) [arXiv:hep-th/0611304]. 

[11] A. Ashtekar and B. Krishnan, "Isolated and dynamical horizons and their applications," 
Living Rev. Rel. 7, 10 (2004) [arXiv:gr-qc/0407042] 

[12] S. de Haro, S. N. Solodukhin and K. Skenderis, "Holographic reconstruction of spacetime 
and renormalization in the AdS/CFT correspondence," Commun. Math. Phys. 217, 595 

(2001) [arXiv:hep-th/0002230]. 

K. Skenderis, "Lecture notes on holographic renormalization," Class. Quant. Grav. 19, 5849 

(2002) [arXiv:hep-th/0209067]. 

[13] L.D. Landau and E.M. Lifshitz, Fluid Mechanics, Pergamon Press, NewYork, 1987, 2nd ed. 



12 



